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Abstract 

Meteorological measurements in the Walnut Gulch catchment in Arizona were used to synthesize a distributed, 
hourly-average time series of data across a 26.9 by 12.5 km area with a grid resolution of 480 m for a continuous 
18-month period which included two seasons of monsoonal rainfall. Coupled surface-atmosphere model runs 
established the acceptability (for modelling purposes) of assuming uniformity in all meteorological variables 
other than rainfall. Rainfall was interpolated onto the grid from an array of 82 recording rain gauges. These 
meteorological data were used as forcing variables for an equivalent array of stand-alone Biosphere-Atmosphere 
Transfer Scheme (BATS) models to describe the evolution of soil moisture and surface energy fluxes in response 
to the prevalent, heterogeneous pattern of convective precipitation. The calculated area-average behaviour was 
compared with that given by a single aggregate BATS simulation forced with area-average meteorological data. 
Heterogeneous rainfall gives rise to significant but partly compensating differences in the transpiration and the 
intercepted rainfall components of total evaporation during rain storms. However, the calculated area-average 
surface energy fluxes given by the two simulations in rain-free conditions with strong heterogeneity in soil mois- 
ture were always close to identical, a result which is independent of whether default or site-specific vegetation 
and soil parameters were used. Because the spatial variability in soil moisture throughout the catchment has the 
same order of magnitude as the amount of rain falling in a typical convective storm (commonly 10% of the veg- 
etation’s root zone saturation) in a semi-arid environment, non-linearity in the relationship between transpira- 
tion and the soil moisture available to the vegetation has limited influence on area-average surface fluxes. 


Introduction 

This research note addresses the question, ‘In the 
extreme case of the semi-arid environment subject to 
convective rainfall, are the ensuing real patterns of soil 
moisture such that they have significant impact on the 
area-average description of vegetation in meteorological 
models?’ It is best read as a sequel to Arain et al. (1996). 

Arain et al. investigated the adequacy of simple rules 
for defining the effective, area-average (or aggregate) 
value of the vegetation-related parameters required to 
specify surface-atmosphere interactions for heterogeneous 
mixes of vegetation covers. They used data from the 
First International Satellite Land Surface Climatology 
Field Experiment (FIFE; Sellers et al., 1988) to validate 
and initiate the BATS-ABL coupled model of interac- 


tions between land surfaces and the atmospheric bound- 
ary layer (Arain et al ., 1996). In most conditions they 
found that so-called aggregation rules (Shuttleworth, 
1992) can provide area-average values of vegetation 
parameters which, tvhen applied within the Biosphere- 
Atmosphere Transfer Scheme (BATS; Dickinson et al . , 
1986), calculate values of surface fluxes that are accept- 
ably similar to those calculated with explicit recognition 
of different patches of vegetation cover in the coupled 
model. 

As a side-issue in their paper, Arain et al. (1996) dis- 
covered that when they included a significant proportion 
of irrigated crops amongst heterogeneous vegetation in an 
otherwise dry landscape, the performance of the hypo- 
thetical aggregation rules was compromised. This result 
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is explained by the non-linear response of vegetation to 
changing soil moisture in the BATS and, it is assumed, 
in nature as well — see Arain et ai (1996) for greater 
detail. 

An analysis of satellite-derived land cover classes 
shows that the proportion of irrigated land falling within 
an individual grid square used in meteorological models 
is rarely greater than 10% in the USA. (The USA is the 
only region for which such data sets are readily available 
at this time.) Moreover, coupled modelling studies using 
the BATS-ABL show that including irrigated vegetation 
in this proportion still allows the use of aggregation 
rules with acceptable accuracy in the simulated area- 
average fluxes. Hence, the presence of patches of 
artificially wetted (irrigated) soil is not a serious issue. 

However, Arain et ai (1996) suggested that in certain 
circumstances extreme differences in soil moisture might 
occur naturally and they recommended that a study be 
made to determine if poor aggregation rule performance 
exists in natural situations. This note explores this 
specific issue for the deliberately extreme case of soil 
moisture heterogeneity resulting from intermittent local- 
ized convective rain falling in a semi-arid region of the 
southwestern USA. 

Field Data 

The data used in this study are from the Walnut Gulch 
Experimental Watershed (31° 43' N, 110° 00' W), which 
is located 120 km southeast of Tucson, Arizona, and 
which comprises the upper 150 km 2 of the Walnut Gulch 
drainage basin (see Fig. 1). Here, soil types range from 
clays and silts to well-cemented boulder conglomerates 
(Renard et ai , 1993), with gravelly and sandy loam sur- 


face textures containing an average of 30% rock and little 
organic matter (Kustas and Goodrich, 1994). The topog- 
raphy comprises gently rolling hills incised by rather 
steep drainage channels. The mixed grass-brush range- 
land vegetation, which is typical of southeastern Arizona 
and southwestern New Mexico, ranges from 20 to 60% 
in coverage. 

The catchment receives 250-500 mm of precipitation 
annually, two-thirds of it falling as convective precipita- 
tion during a summer monsoon season encompassing the 
months of July through to mid-September. The balance 
of precipitation falls mainly during winter frontal storms 
of Pacific origin. The monsoonal precipitation is usually 
from high-intensity convective thunderstorms of limited 
areal extent. The average run-off in the ephemeral 
streams is small (~5%) and of short duration, which 
implies that one-dimensional models, like the BATS, can 
provide a reasonably realistic description of this region's 
surface-atmosphere interactions. The 85 recording rain 
gauges which are currently maintained in and around the 
catchment — an average of one rain gauge per 1.76 km 2 , 
along with 11 runoff measuring flumes, and two sets of 
weather and energy flux measuring stations are indicated 
in Fig. 1. Rainfall and runoff are recorded continuously 
on analogue charts and digitized into breakpoint data. 
Practical temporal resolution of the rainfall and runoff 
data is ± 5 min, and ± 2 min., respectively (Renard et 
a/., 1993; Kustas and Goodrich, 1994). Weather and 
energy flux measuring stations are recorded on data log- 
gers every 20 minutes (Stannard et ai , 1994; Kustas and 
Goodrich, 1994). For additional information on the Wal- 
nut Gulch catchment and the field systems deploved 
there see Kustas et ai (1991), (1994a), (1994b), (1994c), 
Kustas and Goodrich (1994), Stannard et ai (1994), 



Fig. 1 The area described by the distributed array of stand-alone BATS models which includes the Walnut Gulch catchment in Arizona 
The figure shows the catchment boundary the location of rain gauges, flumes and the Kendall and Lucky Hills nucrometeorological stations. 
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Moran et ai (1994), Humes et ai (1994), Weltz et ai 
(1994), Menenti and Ritchie (1994) and Hipps et ai 
(1994). 

SYNTHESIS OF DISTRIBUTED DATA ARRAYS 

This modelling study commenced on July 16, 1990 and 
continued through to December 31, 1991. It therefore 
included two seasons of monsoon rainfall, the first season 
being wetter than average and the second season drier 
than average. The modelling strategy (described below) 
requires a two-dimensional, distributed array of hourly- 
average near-surface meteorological forcing variables, for 
the stand-alone version of the BATS for the 18-month 
study period. These data were synthesized from available 
measurements and projected on to a 56 column by 26 
row array that covers the Walnut Gulch catchment at a 
480 m grid resolution (Figure 1). The 651 grid squares, 
which are located within the boundary of the Walnut 
Gulch catchment, are those for which model comparisons 
were made in this study. 

Comparisons between the meteorological observations 
at the two weather and energy flux measuring stations 
showed that, when it was not raining, the inter-site dif- 
ferences in solar radiation, temperature and vapour pres- 
sure deficit were low (typically less than 20 Wm 2 , 1 °C, 
and 0.14 kPa respectively), even when the difference in 
the observed soil moisture at these sites was at a max- 
imum. Larger, short-term differences were observed 
(particularly in temperature and vapour pressure deficit) 
while it was raining, if the intensity and, in particular, 
the timing of rainfall differed between the two sites. 
Large short term differences in solar radiation were typi- 
cal on cloudy days. However, model sensitivity tests, 
using the two observed time series of weather variables 
as alternate forcings, demonstrated that the overall 
results of this study are not sensitive to this meteorologi- 
cal variability. This is because the periods for which 
there were significant differences in the meteorological 
variables between the two sites is a small proportion of 
the time for which model integrations were made. More- 
over, studies with the BATS-ABL coupled model 
(described later) show' that area-average surface fluxes 
have little sensitivity to atmospheric advection that might 
have been generated by the most extreme observed dif- 
ferences in soil moisture. 

In consequence of the just-described observations and 
modelling results, it is appropriate to assume that the 
values of incoming solar radiation, near-surface air tem- 
perature and vapour pressure deficit required by the 
stand-alone version of the BATS are, in effect, constant 
across the 26.9 by 12.5 km study area. For reasons of 
data consistency, the preferred values of meteorological 
forcing variables applied during modelling were those from 
the Kendall micrometeorological station (Fig. 1). How- 
ever, for about 1% of the time, data were missing from 


this site so data from Lucky Hills were substituted, and 
for a few r brief periods (less than 1% of the time) when 
neither station was operational, it was necessary to use 
data from a nearby Fort Huachuca meteorological sta- 
tion. Tests made at times when data from both Walnut 
Gulch and Fort Huachuca were available confirmed that 
differences in these forcings were insignificant for mod- 
elling purposes. 

Distributed precipitation fields were provided by inter- 
polation from the rain gauge network. To ensure consis- 
tency, interpolation was made from the 82 gauges for 
which continuous rainfall data were available for the 
whole study period. Raw' rain gauge data were digitized 
from analogue charts and then converted into an hourly 
time series for the 18-month study period. Previous 
studies with this rain gauge network (Haider, 1994) have 
demonstrated that the multiquadric-biharmonic tech- 
nique is the preferred method of interpolating precipita- 
tion because it combines practicality with acceptable 
accuracy. An hourly time series of precipitation data w^as 
synthesized for each node in the distributed data grid by 
applying this technique at each hour of the study period 
to interpolate the rain gauge data to the equivalent grid 
location. 

BATS MODEL CALIBRATION 

To test the sensitivity of the specific aggregated results 
to the vegetation parameters in the BATS, the aggregate 
modelling tests (described later) w'ere performed using 
both the default BATS semi-arid vegetation parameters 
and the site-specific parameter values observed at the 
Walnut Gulch study site. 

Performance of these two parameter sets was evaluated 
by comparing the resulting ^modelled surface fluxes with 
the measurements at the two observation sites. However, 
the value of this comparison was restricted by known 
weaknesses in the flux measuring instruments used at 
these sites. The observation of sensible heat flux was 
made at a height of 9 m, using the eddy correlation 
method with a single l GilF propeller to measure the ver- 
tical wind velocity (Stannard et al , 1994). A previous 
comparison w r ith relevant energy- and water-balance data 
(Williams, 1996) indicates that the routine use of this 
wind sensor results, on average, in about a 30% under- 
measurement of sensible heat flux. Therefore, a fixed 
correction of this size was applied to the measured sensi- 
ble heat flux prior to its use in the energy budget, to cal- 
culate latent heat flux for the purpose of comparison 
with model calculations. 

The standard BATS parameters appropriate to the 
Walnut Gulch study site corresponded to surface 
cover/vegetation class number 11, the semi-desert 
default (from the 18 available classes), soil textual class 3 
(from 12 textural classes ranging from sand [1] to clay 
[12]), and soil colour class 2 (from the 8 available colour 
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classes ranging from light [1] to dark [8]) (Dickinson et 
a!., 1986). Simulations using these parameters provided 
poor agreement with observed surface energy fluxes, both 
in terms of the overall magnitude of the fluxes (latent 
heat fluxes being on average underestimated by 58%, 
while sensible heat fluxes were overestimated by 33%), 
and in terms of the modelled daily pattern (the model 
calculated that latent heat fluxes fell rapidly after an early 
morning peak, while observations suggested a more uni- 
form daily pattern). 

Selection of site-specific parameters was largely based 
on the literature, much of w%ch resulted from the Mon- 
soon-90 field campaign which was held in the Walnut 
Gulch catchment from July 23 to August 10, 1990, and 
which has been documented extensively (e.g. Kustas et 
al. y 1991, 1994a, 1994b, 1994c, Kustas and Goodrich, 
1994, Stannard et al . , 1994, Moran et al ., 1994, Humes et 
ul •, 1994, Weltz et al . , 1994, Menenti and Ritchie, 1994, 
Hipps et al . , 1994). These sources allowed the confident 
specification of many BATS vegetation parameters 
(including percentage vegetation cover, roughness length, 
zero-plane displacement height, albedo, and maximum 
and minimum leaf area indices), and many soil parame- 
ters (including porosity, saturated hydraulic conductivity, 
and wet and dry soil albedos). In addition, a prior mod- 
elling study carried out at a Sonoran Desert site near 
Tucson, Arizona (Unland et al., 1996) suggested that two 
other parameters (the minimum stomatal resistance of 
the vegetation and the fraction of water extracted from 
the upper soil layer) require revision if the BATS is to 
provide realistic simulation of semi-arid vegetation. The 
values of these two parameters were therefore deter- 
mined by a simple optimization process, against the 
observed surface energy fluxes in an analogous way to 
Unland et ai (1996). 

Use of the revised, site-specific parameters greatly 
enhanced the quality of the comparison between the 
BATS modelled fluxes and observed fluxes, resulting in 
correlation coefficients of 0.875, 0.912, and 0.970 and 
root mean squared errors of 41.7 Wm -2 , 50.8 Wm' 2 , 
and 20.3 Wm 2 for latent, sensible and ground heat flux 
respectively. 

Modelling Studies 

This project explored two modelling issues. The first was 
to check the validity of assuming uniform meteorology 
(apart from precipitation) across the Walnut Gulch 
catchment, while the second appraised the impact of het- 
erogeneous precipitation on vegetation-atmosphere inter- 
actions. 

COUPLED MODELLING TESTS 

The BATS-ABL model was used to investigate the 
significance of possible patch-to-patch atmospheric 
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advection in the presence of extreme soil moisture het- 
erogeneity observed across the study site. An initial run 
w'as made using the distributed array of stand-alone 
BATS models forced with the distributed array of mete- 
orological data (see below) to determine the most 
extreme differences in modelled soil moisture status in 
the catchment. The soil moisture status and soil and veg- 
etation temperatures for the two most extreme elements 
in the model array were selected and used as initial val- 
ues for moist and dry patches in BATS-ABL, while the 
area-average meteorology at the time these extreme val- 
ues occurred was used to initiate the atmospheric vari- 
ables in that model (Arain, 1994), Arain et al . , 1996). 

The BATS-ABL was run from these initial conditions 
with a modelled domain size of 2000 m and with alter- 
nate vegetation patches (with wet and dry soil) each 1000 
m wide (Arain et ai , 1996). On the basis of the surface- 
atmosphere interaction simulated with the BATS-ABL, 
area-average values of the surface exchanges and of the 
near-surface meteorological variables at 2 m were calcu- 
lated across the modelled domain. A separate calculation 
of the area-average surface exchanges was then made, 
with the area-average value of the meteorological vari- 
ables at 2 m used as forcing variables (rather than the 
spatially varying values of forcing variables that had been 
modelled to occur at this height in the BATS-ABL run). 
In this way, an assessment could be made of the accept- 
ability of assuming a uniform meteorology (in all but 
precipitation) across the Walnut Gulch catchment in the 
distributed array model runs. The results of this coupled 
modelling study were reassuring in that the calculated 
area-average fluxes, given by the BATS-ABL (which 
acknowledge spatial variability in near-surface meteorol- 
ogy coupled to variability in soil moisture), agreed with 
those calculated using an enforced uniformity in near 
surface meteorology to within 11.5%, 8.6% and 7.1% for 
latent, sensible and ground heat flux respectively in these 
conditions of most extreme soil moisture difference. 

AGGREGATION MODELLING STUDIES 

Two stand-alone BATS modelling studies were used to 
assess the impact of heterogeneous precipitation on 
vegetation-atmosphere interactions. An individual (aggre- 
gate) BATS model and a distributed array of single point 
BATS models were created and parameterized identically 
over the Walnut Gulch catchment. The distributed 
BATS model was driven with the spatially distributed 
meteorological forcing data derived as described above, 
while the aggregate BATS model was forced with meteo- 
rological data which was the area-average of that used 
for the distributed array for each hour. (Recall that, in 
practice, precipitation is the only variable considered to 
have significant spatial heterogeneity). 

The calculated surface exchanges for the array of 
single point BATS models were averaged and compared 
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Fig. 2 Comparison between the hourly-average (a) the latent-heat flux ; (b) sensible-heat flux; and (c) ground heat flux as calculated by 
the single, aggregate BATS model forced with area-average meteorological data and the distributed array of BATS models forced with dis- 
tributed data for all rain-free hours in the 18 month study period 


to the output from the aggregate BATS model which, it 
was hypothesized, provided an adequate area-average 
description. In this way, the area-average interaction 
with the realistic representation of precipitation-induced 
heterogeneity in soil moisture can be compared with that 
given when heterogeneity is neglected in the modelling 
process. If the two outputs disagree with one another, 
then failure of aggregation rules (analogous to that which 
Arain et al. (1996) reported as possible with patches of 
irrigated agriculture) can be deemed to have occurred in 
these albeit extreme, but nonetheless natural, conditions. 

The results indicate no evidence of significant aggrega- 
tion rule failure with naturally-occurring differences in 
soil moisture. Modelling results are presented here only 
for runs using the site-specific BATS parameters but 
simulations were also made with the BATS default 
parameters (see above). Although the area-average sur- 
face energy fluxes, calculated with the default parameter- 
ization in both the aggregate and the distributed array 
BATS models, differed from the equivalent runs using 
the site-specific parameters, none of the conclusions 
given below is affected because these conclusions are 
based on the difference between aggregate and dis- 
tributed array runs made with the same set of model 
parameters. 

Figure 2 shows a comparison between the calculated 
latent heat, sensible heat and ground heat fluxes given by 
the aggregate BATS model and those given by the dis- 
tributed array of BATS models for hours when there 
was no intercepted water on the vegetation canopy dur- 
ing the 18 month study. Clearly surface energy fluxes 
agree extremely well, to within a few 10s Wm 2 — this 
being the criterion for success adopted by Arain et al 
(1996). 

Notwithstanding this success in rain-free conditions, 
an interesting side-effect of spatial variability in precipi- 
tation is apparent in this study, and is related to the 
interception of rainfall by the vegetation canopy. When 


the rain falling in a storm covers a large proportion 
(approximately 60% or more) of the study area modelled 
there is an inconsequential difference between the calcu- 
lated components of evapouration for simulations 
made with the aggregate BATS model and the dis- 


Transplration 



DOY 1990 (DOY 213 is August 1,1990) 


Evaporation of Intarcepted Water 



DOY 1990 (DOY 213 is August 1,1990) 


Interpolated Precipitation 



DOY 1990 (DOY 213 is August 1,1990) 


Fig. 3 Time series of calculated rates of (a) transpiration , (b) 
evaporation of intercepted water as given by the single, aggregate 
BATS model forced with area-average meteorological data (full 
line) and the distributed array of BATS models forced with dis- 
tributed data (broken line) and (c) the interpolated precipitation 
from (DOY) day of year 213 to 215 of 1990 (August 1 -August 3, 
1990). There were 7 localized rain storms during this period. 
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Table I. Modelled Water Balance Components for the Aggre- 
gate and Distributed BATS models using site-specific parame- 
ters. 


Modelled Water Balance from July 16, 1990-December 31, 
1991 (534 Day Model Run) 


Cumulative water 
balance component 

Aggregate 

model 

(mm) 

Distributed 

model 

(mm) 

% Difference* 

Evaporation flux 

446.27 

467.29 

0.22% 

Precipitation 

426.46 

426.34 

-0.03% 

Runoff 

21.33 

20.33 

-4.92% 

A Storage from soil 

-54.53 

-55.05 

0.99% 

moisture 




Transpiration 

359.41 

366,05 

1.81% 

Interception 

25.91 

19.51 

-32.80% 

Soil evaporation 

73.62 

74.59 

1.30% 


W ater balance closure 



X Input 

1 Output 

X in/X Out 


(mm) 

(mm) 


Aggregate 

459.66 

466.27 

98.6% 

Distributed 

461.06 

467.29 

98.7% 


* [(Distributed - Aggregate)/Distributed] X 100% 


tributed array of BATS model. However, when storms 
of limited extent occur, this is not the case. Figure 3 
shows an example of the area-average, model-calculated 
evaporation flux for the study area that originate (a) as 
transpiration, and (b) as the evaporation of intercepted 
water, along with (c) the interpolated rainfall for 3 days 
when there were 7 short storms, 3 of which covered 
approximately 25% of the modelled area whereas the 
remaining 4 covered less than 10% of the modelled area. 
Inspection of Fig. 3 reveals that localized storms of this 
type give rise to differences in the modelled area-average 
rate of evaporation of intercepted water between the 
aggregate and distributed array models, which can be as 
large as 30-50% for short periods. 

This difference arises because the aggregate model 
rainfall is (wrongly) assumed to be equally distributed 
across the modelled area, so all elements in the grid con- 
tribute to the net interception loss, and its modelled 
interception rate is therefore higher. At the same time, 
there is (in part) a compensating difference in the 
area-average transpiration flux, because it is suppressed 
for the whole study area and not just for the rained-on 
portion. However, the compensation between these two 
components of evaporation loss is not perfect and, in 
general, the erroneous enhancement of calculated inter- 
ception loss for the aggregate BATS model slightly 
exceeds the erroneous reduction in the calculated tran- 
spiration loss during the storm. Consequently, more 
water evaporates and less rain infiltrates during storms. 


causing the modelled soil moisture and transpiration to 
be slightly less on rain-free days. This small reduction in 
inter-storm, area-average transpiration is 1.8% as mod- 
elled. Table 1 lists this difference, as well as the differ- 
ences between the cumulative water balance components 
for both simulations. 

Bare soil evaporation rate is known significant^ to 
enhance evaporation rates for sparse, semi-arid crops for 
a few days after each rain storm (Gash et al, 1991). The 
presence of this enhanced soil evaporation, which is 
restricted to the rained-on portion of the landscape, 
could generate differences between calculations made 
with the distributed and aggregate models for reasons 
which are essentially identical to those for the intercep- 
tion process described above. Such short-term differ- 
ences can occur for storms which cover only part of the 
modelled area, but the cumulative effect over the whole 
study period of differences in soil evaporation calculated 
with the two models is small (Table 1). 

Discussion and Conclusions 

On the basis of the above comments, it is clear that het- 
erogeneous rainfall can result in a poor description of 
short-term evaporation rates during and immediately 
after rainfall because the description of the rainfall inter- 
ception and soil evaporation processes is imperfect in a 
one-dimensional, aggregate representation. This feature 
has been recognized for some time, (Shuttleworth, 1988). 
However, in the case of semi-arid regions, where there is 
little runoff and where water availability is the limiting 
control on evaporation, maintaining the water balance 
necessarily requires that modelled transpiration is 
reduced between rain storms to compensate for this addi- 
tional within-storm loss due to the evaporation of inter- 
cepted rain. In the case of soil evaporation, there is 
similar compensation for patches of bare soil, between 
greater soil evaporation calculated with the distributed 
model immediately after rainfall and reduced soil evapo- 
ration calculated between storms. Such compensation 
between the separate components of total evaporation 
will not necessarily occur in regions where evaporation 
rates are energy limited and runoff is plentiful. 

Notwithstanding this last mentioned complication, the 
primary focus of the present study was to investigate 
whether realistic patterns of soil moisture have significant 
impact on the area-average description of vegetation in 
meteorological models in the extreme case of a semi-arid 
environment subject to intermittent but intense convec- 
tive rainfall. The results are definitive in this respect. In 
rain-free conditions but with natural, rain-induced, het- 
erogeneous soil moisture, the difference between the cal- 
culated surface energy fluxes given by a single, aggregate 
BATS model and that given by a distributed arrav of 
BATS models is very small indeed. 

The basic reason for this is that the amount of rain in 
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an individual precipitation event is always a fairly small 
proportion of the total amount of water that can be 
stored within the rooting zone of the simulated vegeta- 
tion. In this study, the area-average soil moisture stored 
in the rooting zone varies with season: typically, the soil 
moisture store is about half full during the monsoon sea- 
son, and close to empty after an extended dry period. 
However, because local convective storms occur ran- 
domly in space, spatial variations in available soil mois- 
ture across the landscape at any time have the same 
order of magnitude as the amount of rain in one typical 
rainfall event. Convective rain storms can be heavy and 
localized in this environment — the largest storm observed 
at any grid point in the modelled array during the study 
period was about 50 mm, for instance, but the amount of 
water that can be stored in the rooting zone of the vege- 
tation (336 mm) is large in comparison with this. Hence, 
at any time, only about 10% of the non-linear relation- 
ship between soil moisture and evaporation is being sam- 
pled within a given soil moisture pattern, so that the 
non-linear relationship can be considered approximately 
linear over this small, sampled range. Given that this is 
the case in this extreme semi-arid situation, it is likely 
that it is more generally true elsewhere. 
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